Workflow 3: comparing observations to emissions

In addition to observation files, ancillary data can also be added to an openghg object store which can be used to perform analysis.

At the moment, the accepted files include:

• Footprints - regional outputs from an LPDM model (e.g. NAME)

• Emissions/Flux - estimates of species emissions within a region

• [+Boundary conditions - to be added+]

• Global CTM output (e.g. GEOSChem)

These inputs must adhere to an expected format and are expected to minimally contain a fixed set of inputs.

At the moment, the expected format for these files is created through standard methods from within the ACRG repository.

Loading data sources into the object store

For this tutorial we will again set up a temporary object store to store our data.

See 1_Adding_observation_data.ipynb for more details and advice on how to create a more permanent object store. Once a permanent object store is set up, these steps would only need to be performed once. Any added data can then be retrieved using searches.

import os
import tempfile

tmp_dir = tempfile.TemporaryDirectory()
os.environ["OPENGHG_PATH"] = tmp_dir.name   # temporary directory

For this, we will add observation, footprint and flux data to the object store. This data relates to Tacolneston (TAC) site within the DECC network and the area around Europe (EUROPE domain).

from openghg.client import process_obs

obs_file = "../data/DECC/tac.picarro.1minute.100m.201208.dat"

site="tac" 
network="DECC"
height="100m"
species="ch4"

process_obs(files=obs_file, data_type="CRDS", site=site, network=network, inlet=height)

  0%|          | 0/1 [00:00<?, ?it/s]

Processing: tac.picarro.1minute.100m.201208.dat:   0%|          | 0/1 [00:00<?, ?it/s]

Processing: tac.picarro.1minute.100m.201208.dat: 100%|██████████| 1/1 [00:00<00:00,  2.75it/s]

Processing: tac.picarro.1minute.100m.201208.dat: 100%|██████████| 1/1 [00:00<00:00,  2.75it/s]

{'processed': {'tac.picarro.1minute.100m.201208.dat': {'ch4': '4a5d62d5-4065-49bf-aad8-c8c7b86e72e8',
   'co2': 'dc7a36ce-064b-4e2d-8464-8e3151cbdab5'}}}

from openghg.client import process_footprint

fp_file_path = "../data/footprints/TAC-100magl_EUROPE_201208.nc"

site="tac" 
height="100m"
domain="EUROPE"
model="NAME"

process_footprint(files=fp_file_path, site=site, height=height, domain=domain, model=model)

{'tac_europe_NAME_100m': 'a571d49f-f769-4d36-bc7c-088e6ac15c70'}

from openghg.util import retrieve_example_data

flux_data = retrieve_example_data(path="flux/ch4-ukghg-all_EUROPE_2012.tar.gz")

from openghg.util import retrieve_example_data
from openghg.client import process_flux

flux_data = retrieve_example_data(path="flux/ch4-ukghg-all_EUROPE_2012.tar.gz")

site="tac" 
domain="EUROPE"
date = "2012"

source_waste = "waste"
source_energyprod = "energyprod"

flux_data_waste = [filename for filename in flux_data if source_waste in str(filename)][0]
flux_data_energyprod = [filename for filename in flux_data if source_energyprod in str(filename)][0]

process_flux(files=flux_data_waste, species=species, source=source_waste, domain=domain, date=date)
process_flux(files=flux_data_energyprod, species=species, source=source_energyprod, domain=domain, date=date)

{'ch4_energyprod_europe_2012': 'bbb3dc2e-4cfc-45f2-a26e-ffb41009a4ba'}

Creating a model scenario

With this ancillary data, we can start to make comparisons between model data, such as bottom-up inventories, and our observations. This analysis is based around a ModelScenario class which can be created to link together observation, footprint and emissions data.

Boundary conditions and other model data will be added soon

Above we loaded observation data from the Tacolneston site into the object store. We also added an associated footprint (sensitivity map) and anthropogenic emissions maps both for a domain defined over Europe.

To access and link this data we can set up our ModelScenario instance using a similiar set of keywords. In this case we have also limited ourselves to a date range:

from openghg.analyse import ModelScenario

start_date = "2012-08-01"
end_date = "2012-09-01"

scenario = ModelScenario(site=site, 
                         inlet=height, 
                         domain=domain, 
                         species=species, 
                         source=source_waste, 
                         start_date=start_date,
                         end_date=end_date)

/opt/hostedtoolcache/Python/3.9.10/x64/lib/python3.9/site-packages/xarray/core/indexes.py:97: FutureWarning: Indexing a timezone-naive DatetimeIndex with a timezone-aware datetime is deprecated and will raise KeyError in a future version.  Use a timezone-naive object instead.
  indexer = index.slice_indexer(

Adding obs_surface to model scenario
Updating any inputs based on observation data
site: tac, species: ch4, inlet: 100m

Adding footprint to model scenario
Adding flux to model scenario

Using these keywords, this will search the object store and attempt to collect and attach observation, footprint and flux data. This collected data will be attached to your created ModelScenario. For the observations this will be stored as the ModelScenario.obs attribute. This will be an ObsData object which contains metadata and data for your observations:

scenario.obs

ObsData(data=<xarray.Dataset>
Dimensions:                    (time: 20509)
Coordinates:
  * time                       (time) datetime64[ns] 2012-08-01T00:00:30 ... ...
Data variables:
    mf                         (time) float64 1.915e+03 1.912e+03 ... 1.942e+03
    mf_variability             (time) float64 1.431 2.018 1.767 ... 0.865 0.74
    mf_number_of_observations  (time) float64 20.0 19.0 20.0 ... 19.0 20.0 6.0
Attributes: (12/24)
    data_owner:            Simon O'Doherty
    data_owner_email:      s.odoherty@bristol.ac.uk
    inlet_height_magl:     100m
    comment:               Cavity ring-down measurements. Output from GCWerks
    long_name:             tacolneston
    site:                  tac
    ...                    ...
    sampling_period_unit:  s
    station_longitude:     1.13872
    station_latitude:      52.51775
    station_long_name:     Tacolneston Tower, UK
    station_height_masl:   50.0
    scale:                 WMO-X2004A, metadata={'site': 'tac', 'instrument': 'picarro', 'sampling_period': '60', 'inlet': '100m', 'port': '9', 'type': 'air', 'network': 'decc', 'species': 'ch4', 'calibration_scale': 'wmo-x2004a', 'long_name': 'tacolneston', 'data_owner': "Simon O'Doherty", 'data_owner_email': 's.odoherty@bristol.ac.uk', 'station_longitude': 1.13872, 'station_latitude': 52.51775, 'station_long_name': 'Tacolneston Tower, UK', 'station_height_masl': 50.0, 'inlet_height_magl': '100m', 'data_type': 'timeseries', 'comment': 'Cavity ring-down measurements. Output from GCWerks', 'conditions_of_use': 'Ensure that you contact the data owner at the outset of your project.', 'source': 'In situ measurements of air', 'Conventions': 'CF-1.8', 'file_created': '2022-03-22 10:42:54.957441+00:00', 'processed_by': 'OpenGHG_Cloud', 'sampling_period_unit': 's', 'scale': 'WMO-X2004A'})

To access the undelying xarray Dataset containing the observation data use ModelScenario.obs.data:

ds = scenario.obs.data

The ModelScenario.footprint attribute contains the linked FootprintData (again, use .data to extract xarray Dataset):

scenario.footprint

FootprintData(data=<xarray.Dataset>
Dimensions:               (time: 372, lon: 391, lat: 293, lev: 1, height: 20)
Coordinates:
  * time                  (time) datetime64[ns] 2012-08-01 ... 2012-08-31T22:...
  * lon                   (lon) float32 -97.9 -97.55 -97.2 ... 38.68 39.03 39.38
  * lat                   (lat) float32 10.73 10.96 11.2 ... 78.59 78.82 79.06
  * lev                   (lev) object 'From     0 -    40m agl'
  * height                (height) float32 500.0 1.5e+03 ... 1.85e+04 1.95e+04
Data variables:
    fp                    (lat, lon, time) float32 0.0 0.0 ... 7.629e-06
    temperature           (time) float32 15.64 14.35 14.01 ... 14.08 13.07 13.31
    pressure              (time) float32 994.9 994.3 ... 1.01e+03 1.009e+03
    wind_speed            (time) float32 8.483 8.844 8.768 ... 2.48 3.621 4.743
    wind_direction        (time) float32 193.2 203.1 201.4 ... 42.42 117.3 172.3
    PBLH                  (time) float32 510.2 341.6 291.5 ... 132.5 92.75 92.75
    release_lon           (time) float32 1.139 1.139 1.139 ... 1.139 1.139 1.139
    release_lat           (time) float32 52.52 52.52 52.52 ... 52.52 52.52 52.52
    particle_locations_n  (height, lon, time) float32 0.0 3.248e-05 ... 0.0 0.0
    particle_locations_e  (height, lat, time) float32 0.0 0.0 0.0 ... 0.0 0.0
    particle_locations_s  (height, lon, time) float32 0.0 0.0 0.0 ... 0.0 0.0
    particle_locations_w  (height, lat, time) float32 0.0 0.0 0.0 ... 0.0 0.0
Attributes: (12/16)
    author:           OpenGHG Cloud
    processed:        2022-03-22 10:42:55.132208+00:00
    data_type:        footprints
    site:             tac
    height:           100m
    domain:           europe
    ...               ...
    min_longitude:    -97.9
    max_latitude:     79.057
    min_latitude:     10.729
    time_resolution:  standard_time_resolution
    heights:          [  500.  1500.  2500.  3500.  4500.  5500.  6500.  7500...
    variables:        ['fp', 'temperature', 'pressure', 'wind_speed', 'wind_d..., metadata={'data_type': 'footprints', 'site': 'tac', 'height': '100m', 'domain': 'europe', 'model': 'name', 'start_date': '2012-08-01 00:00:00+00:00', 'end_date': '2012-08-31 22:00:00+00:00', 'max_longitude': 39.38, 'min_longitude': -97.9, 'max_latitude': 79.057, 'min_latitude': 10.729, 'time_resolution': 'standard_time_resolution', 'heights': [500.0, 1500.0, 2500.0, 3500.0, 4500.0, 5500.0, 6500.0, 7500.0, 8500.0, 9500.0, 10500.0, 11500.0, 12500.0, 13500.0, 14500.0, 15500.0, 16500.0, 17500.0, 18500.0, 19500.0], 'variables': ['fp', 'temperature', 'pressure', 'wind_speed', 'wind_direction', 'pblh', 'release_lon', 'release_lat', 'particle_locations_n', 'particle_locations_e', 'particle_locations_s', 'particle_locations_w']}, flux={}, bc={}, species='NA', scales='FIXME', units='FIXME')

And the ModelScenario.fluxes attribute can be used to access the FluxData. Note that for ModelScenario.fluxes this can contain multiple flux sources and so this is stored as a dictionary linked to the source name:

scenario.fluxes

{'waste': FluxData(data=<xarray.Dataset>
 Dimensions:  (time: 1, lat: 293, lon: 391)
 Coordinates:
   * time     (time) datetime64[ns] 2012-01-01
   * lat      (lat) float32 10.73 10.96 11.2 11.43 ... 78.36 78.59 78.82 79.06
   * lon      (lon) float32 -97.9 -97.55 -97.2 -96.84 ... 38.32 38.68 39.03 39.38
 Data variables:
     flux     (time, lat, lon) float32 0.0 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
 Attributes: (12/17)
     Processed by:     cv18710@bristol.ac.uk
     Processed on:     2020-11-06 11:55:50.312453+00:00
     Raw file used:    /home/cv18710/work_shared/Gridded_fluxes/CH4/UKGHG/uk_f...
     author:           OpenGHG Cloud
     species:          ch4
     domain:           europe
     ...               ...
     max_longitude:    39.38
     min_longitude:    -97.9
     max_latitude:     79.057
     min_latitude:     10.729
     time_resolution:  standard
     frequency:        annual, metadata={'processed by': 'cv18710@bristol.ac.uk', 'processed on': '2020-11-06 11:55:50.312453+00:00', 'raw file used': '/home/cv18710/work_shared/gridded_fluxes/ch4/ukghg/uk_flux_waste_ch4_lonlat_0.01km_2012.nc', 'species': 'ch4', 'domain': 'europe', 'source': 'waste', 'date': '2012', 'author': 'openghg cloud', 'processed': '2022-03-22 10:43:00.404396+00:00', 'start_date': '2012-01-01 00:00:00+00:00', 'end_date': '2012-12-31 23:59:59+00:00', 'max_longitude': 39.38, 'min_longitude': -97.9, 'max_latitude': 79.057, 'min_latitude': 10.729, 'time_resolution': 'standard', 'frequency': 'annual', 'data_type': 'emissions'}, flux={}, bc={}, species='ch4', scales='FIXME', units='FIXME')}

An interactive plot for the linked observation data can be plotted using the ModelScenario.plot_timeseries() method:

scenario.plot_timeseries()

You can also set up your own searches and add this data directly.

from openghg.retrieve import get_obs_surface, get_footprint, get_flux

# Extract obs results from object store
obs_results = get_obs_surface(site=site,
                              species=species,
                              inlet=height,
                              start_date="2012-08-01",
                              end_date="2012-09-01")

# Extract footprint results from object store
footprint_results = get_footprint(site=site,
                                  domain=domain,
                                  height=height,
                                  start_date="2012-08-01",
                                  end_date="2012-09-01")

# Extract flux results from object store
flux_results = get_flux(species=species,
                        domain=domain,
                        source=source_waste)

/opt/hostedtoolcache/Python/3.9.10/x64/lib/python3.9/site-packages/xarray/core/indexes.py:97: FutureWarning:

Indexing a timezone-naive DatetimeIndex with a timezone-aware datetime is deprecated and will raise KeyError in a future version.  Use a timezone-naive object instead.

scenario_direct = ModelScenario(obs=obs_results, footprint=footprint_results, flux=flux_results)

Updating any inputs based on observation data
site: tac, species: ch4, inlet: 100m

You can create your own input objects directly and add these in the same way. This allows you to bypass the object store for experimental examples. At the moment these inputs need to be ObsData, FootprintData or FluxData objects (can be created using classes from openghg.dataobjects) but simpler inputs will be made available.

One benefit of this interface is to reduce searching the database if the same data needs to be used for multiple different scenarios.

Comparing data sources

Once your ModelScenario has been created you can then start to use the linked data to compare outputs. For example we may want to calculate modelled observations at our site based on our linkec footprint and emissions data:

modelled_observations = scenario.calc_modelled_obs()

Caching calculated data

This could then be plotted directlt using the xarray plotting methods:

modelled_observations.plot()  # Can plot using xarray plotting methods

[<matplotlib.lines.Line2D at 0x7f09b4ee4e80>]

To compare the these modelled observations to the ovbservations themselves, the ModelScenario.plot_comparison() method can be used:

scenario.plot_comparison()

The ModelScenario.footprints_data_merge() method can also be used to created a combined output, with all aligned data stored directly within an xarray.Dataset:

combined_dataset = scenario.footprints_data_merge()
combined_dataset

WARNING: calc_bc keyword has not been used to create output
 BoundaryConditions have not been implemented yet.

<xarray.Dataset>
Dimensions:                    (time: 372, lon: 391, lat: 293, lev: 1,
                                height: 20)
Coordinates:
  * time                       (time) datetime64[ns] 2012-08-01T00:00:30 ... ...
  * lon                        (lon) float32 -97.9 -97.55 -97.2 ... 39.03 39.38
  * lat                        (lat) float32 10.73 10.96 11.2 ... 78.82 79.06
  * lev                        (lev) object 'From     0 -    40m agl'
  * height                     (height) float32 500.0 1.5e+03 ... 1.95e+04
Data variables: (12/16)
    mf                         (time) float64 1.922e+03 1.934e+03 ... 1.921e+03
    mf_variability             (time) float64 3.155 2.061 3.597 ... 0.237 0.5685
    mf_number_of_observations  (time) float64 19.56 19.86 19.36 ... 19.68 19.43
    fp                         (lat, lon, time) float32 0.0 0.0 ... 7.629e+03
    temperature                (time) float32 15.64 14.35 14.01 ... 13.07 13.31
    pressure                   (time) float32 994.9 994.3 ... 1.01e+03 1.009e+03
    ...                         ...
    release_lat                (time) float32 52.52 52.52 52.52 ... 52.52 52.52
    particle_locations_n       (height, lon, time) float32 0.0 3.248e-05 ... 0.0
    particle_locations_e       (height, lat, time) float32 0.0 0.0 ... 0.0 0.0
    particle_locations_s       (height, lon, time) float32 0.0 0.0 ... 0.0 0.0
    particle_locations_w       (height, lat, time) float32 0.0 0.0 ... 0.0 0.0
    mf_mod                     (time) float32 37.35 35.34 27.8 ... 14.56 24.52
Attributes: (12/40)
    author:                OpenGHG Cloud
    processed:             2022-03-22 10:42:55.132208+00:00
    data_type:             footprints
    site:                  tac
    height:                100m
    domain:                europe
    ...                    ...
    station_longitude:     1.13872
    station_latitude:      52.51775
    station_long_name:     Tacolneston Tower, UK
    station_height_masl:   50.0
    scale:                 WMO-X2004A
    resample_to:           coarsest

xarray.Dataset

• Dimensions:
  • time: 372
  • lon: 391
  • lat: 293
  • lev: 1
  • height: 20
• Coordinates: (5)
  • time
    (time)
    datetime64[ns]
    2012-08-01T00:00:30 ... 2012-08-...

    array(['2012-08-01T00:00:30.000000000', '2012-08-01T02:00:30.000000000',
           '2012-08-01T04:00:30.000000000', ..., '2012-08-31T18:00:30.000000000',
           '2012-08-31T20:00:30.000000000', '2012-08-31T22:00:30.000000000'],
          dtype='datetime64[ns]')

  • lon
    (lon)
    float32
    -97.9 -97.55 -97.2 ... 39.03 39.38

    units :

    Degrees_east

    array([-97.9  , -97.548, -97.196, ...,  38.676,  39.028,  39.38 ],
          dtype=float32)

  • lat
    (lat)
    float32
    10.73 10.96 11.2 ... 78.82 79.06

    units :

    Degrees_north

    array([10.729, 10.963, 11.197, ..., 78.589, 78.823, 79.057], dtype=float32)

  • lev
    (lev)
    object
    'From 0 - 40m agl'

    array(['From     0 -    40m agl'], dtype=object)

  • height
    (height)
    float32
    500.0 1.5e+03 ... 1.85e+04 1.95e+04

    array([  500.,  1500.,  2500.,  3500.,  4500.,  5500.,  6500.,  7500.,  8500.,
            9500., 10500., 11500., 12500., 13500., 14500., 15500., 16500., 17500.,
           18500., 19500.], dtype=float32)

• Data variables: (16)
  • mf
    (time)
    float64
    1.922e+03 1.934e+03 ... 1.921e+03

    array([1921.95491228, 1934.17964286, 1904.42586207, 1910.21589286,
           1910.98303571, 1916.06807018, 1945.88464286, 1939.17232143,
           1896.24473684, 1923.26982143, 1920.40785714, 1898.48232143,
           1899.695     , 1897.62446429, 1910.28446429, 1920.54982456,
           1904.68482143, 1889.3666    , 1877.90877193, 1879.51310345,
           1881.31053571, 1879.68375   , 1889.87517857, 1901.71910714,
           1915.44034483, 1912.19631579, 1907.82      , 1902.34162162,
           1887.07280702, 1878.59053571, 1882.00803571, 1891.16      ,
           1884.70053571, 1880.84339286, 1898.99631579, 1895.45913793,
           1919.20642857, 1925.26763636, 1923.61285714, 1911.13982759,
           1895.06473684, 1894.23517857, 1881.52267857, 1885.56614035,
           1889.43148148, 1899.71561404, 1892.96321429, 1909.51263158,
           1931.56      , 1931.27807018, 1945.51321429, 1953.29105263,
           1908.47508772, 1897.18245614, 1891.6577193 , 1890.22821429,
           1891.4375    , 1899.59649123, 1888.40722222, 1890.21263158,
           1928.91526316, 1902.87964286, 1908.48719298, 1934.83754386,
           1915.8654386 , 1937.43894737, 1927.93125   , 1925.68375   ,
           1923.60892857, 1925.21947368, 1929.82625   , 1939.80321429,
           1964.37245614, 1970.25333333, 1962.83649123, 1966.64125   ,
           1919.77736842, 1919.33054545, 1895.22755102, 1890.56285714,
    ...
           1892.67482759, 1885.12754386, 1889.70789474, 1886.11611111,
           1888.72964286, 1882.74089286, 1898.81303571, 1930.39789474,
           1948.32982456, 1922.53421053, 1909.30517857, 1913.25877193,
           1901.01875   , 1887.96818182, 1885.98125   , 1887.29298246,
           1898.58892857, 1907.94192982, 1892.74446429, 1890.86714286,
           1908.38741379, 1888.75192982, 1914.90701754, 1977.02438596,
           1899.25693878, 1894.21339286, 1918.315     , 1903.75553571,
           1901.13210526, 1918.88660714, 1922.23767857, 1926.96534483,
           1927.77571429, 1914.22982143, 1908.71540541, 1912.55280702,
           1908.17267857, 1904.51964912, 1901.67821429, 1895.18701754,
           1899.12732143, 1908.02175439, 1906.43789474, 1914.48785714,
           1930.31709091, 1921.56526316, 1923.03350877, 1930.57      ,
           1924.25578947, 1896.51842105, 1895.88053571, 1882.27375   ,
           1881.36140351, 1896.67517857, 1901.145     , 1908.91857143,
           1908.5837931 , 1909.26333333, 1917.44196429, 1914.63660714,
           1899.56071429, 1906.11910714, 1904.53285714, 1883.36614035,
           1884.06444444, 1889.36446429, 1888.92392857, 1887.07839286,
           1886.71526316, 1895.44017857, 1905.53196429, 1911.67160714,
           1901.67196429, 1887.00596491, 1885.28857143, 1888.92267857,
           1893.89196429, 1894.4822807 , 1891.1577193 , 1920.51655172])

  • mf_variability
    (time)
    float64
    3.155 2.061 3.597 ... 0.237 0.5685

    array([ 3.15547368,  2.06126786,  3.59748276,  3.66892857,  0.68776786,
            1.1447193 ,  1.0875    ,  2.42598214,  0.70363158,  0.38753571,
            1.59292857,  0.76003571,  1.79014286,  0.84792857,  0.45126786,
            0.90649123,  1.20183929,  0.60614   ,  0.47289474,  0.45293103,
            0.57214286,  0.93475   ,  1.097     ,  1.16591071,  1.39512069,
            1.66852632,  1.20894643,  1.52262162,  0.90114035,  0.61533929,
            0.46289286,  0.69001786,  0.62707143,  0.72682143,  1.37219298,
            1.02136207,  1.51707143,  1.75078182,  2.19496429,  1.91965517,
            0.61052632,  0.503875  ,  0.88105357,  1.20507018,  0.85846296,
            1.61080702,  0.79258929,  0.75396491,  0.95566   ,  1.17501754,
            4.02826786,  5.11087719,  1.60096491,  0.66722807,  0.55680702,
            0.78816071,  0.41292857,  0.93757895,  0.78711111,  0.36961404,
            1.07963158,  1.14135714,  1.94326316,  4.35326316,  1.64445614,
            1.57629825,  1.30301786,  1.03873214,  0.84682143,  0.76385965,
            1.42810714,  0.843     ,  1.20287719,  3.40861404,  1.33052632,
            1.90723214,  1.65808772,  1.64249091,  0.63722449,  1.07416071,
            0.82524561,  1.88242105,  1.31654386,  0.34094643,  0.347     ,
            0.43610714,  0.56215789,  4.48117857,  0.99878378,  0.74935088,
            0.55873214,  3.24639286,  0.78266667,  1.01933333,  0.45180357,
            0.59094828,  0.50560714,  0.78396429,  0.75445455,  6.14626316,
    ...
            0.31941379,  0.80482456,  0.98540351,  2.82247368,  9.741375  ,
            5.09383929,  0.94755172,  0.42285714,  0.495875  ,  0.72428571,
            0.91493878,  1.48624138,  0.53050877,  1.36371429,  0.64191071,
            1.3465    ,  2.65868966,  0.84281034,  0.72312281,  0.46036842,
            0.62661111,  0.76694643,  1.00001786,  0.93571429,  1.03685965,
            2.08624561,  2.07257895,  0.52169643,  0.62536842,  0.95475   ,
            0.42323636,  0.27360714,  0.24035088,  4.56433929,  1.40229825,
            1.07325   ,  0.95053571,  0.68643103,  0.47635088,  2.72382456,
           14.28482456,  1.187     ,  0.37814286,  0.36403571,  0.39514286,
            0.41545614,  0.58514286,  0.86264286,  0.7557931 ,  1.45476786,
            1.17508929,  1.21256757,  1.88378947,  0.92814286,  0.51735088,
            0.36725   ,  0.35707018,  0.51725   ,  1.13991228,  0.44763158,
            0.40533929,  1.29909091,  1.94638596,  1.71115789,  1.77922807,
            1.36180702,  0.56022807,  0.62519643,  0.98205357,  0.59826316,
            0.56748214,  0.91576   ,  0.91553571,  0.54087931,  0.96394737,
            1.33991071,  1.48116071,  0.93335714,  0.73773214,  0.70771429,
            0.38931579,  0.24991667,  0.28833929,  0.32905357,  0.32771429,
            0.29112281,  0.58935714,  0.44783929,  0.5865    ,  0.5       ,
            0.26677193,  0.40430357,  0.30655357,  0.99175   ,  0.71764912,
            0.23701754,  0.56851724])

  • mf_number_of_observations
    (time)
    float64
    19.56 19.86 19.36 ... 19.68 19.43

    array([19.56140351, 19.85714286, 19.36206897, 19.85714286, 19.875     ,
           19.56140351, 19.75      , 19.51785714, 19.57894737, 19.85714286,
           19.85714286, 19.91071429, 19.89285714, 19.92857143, 19.89285714,
           19.64912281, 19.91071429, 19.62      , 19.64912281, 19.4137931 ,
           19.875     , 19.92857143, 19.94642857, 19.94642857, 19.44827586,
           19.68421053, 19.94642857, 19.56756757, 19.68421053, 19.91071429,
           19.85714286, 19.83928571, 19.91071429, 19.91071429, 19.70175439,
           19.44827586, 19.94642857, 19.96363636, 19.92857143, 19.46551724,
           19.70175439, 19.94642857, 19.98214286, 19.64912281, 19.98148148,
           19.71929825, 19.94642857, 19.71929825, 19.68      , 19.71929825,
           19.96428571, 19.70175439, 19.73684211, 19.70175439, 19.73684211,
           19.96428571, 19.92857143, 19.71929825, 19.97222222, 19.71929825,
           19.71929825, 19.98214286, 19.70175439, 19.71929825, 19.73684211,
           19.70175439, 19.96428571, 19.89285714, 19.64285714, 19.71929825,
           19.98214286, 19.98214286, 19.71929825, 19.71929825, 19.71929825,
           19.96428571, 19.45614035, 19.94545455, 19.97959184, 19.91071429,
           19.73684211, 19.75438596, 19.71929825, 19.98214286, 19.43103448,
           19.92857143, 19.64912281, 19.89285714, 19.51351351, 19.66666667,
           19.875     , 19.89285714, 19.63157895, 19.66666667, 19.92857143,
           19.4137931 , 19.92857143, 19.91071429, 19.89090909, 19.61403509,
    ...
           19.44827586, 19.68421053, 19.68421053, 19.68421053, 19.92857143,
           19.92857143, 19.51724138, 19.94642857, 19.96428571, 20.        ,
           19.97959184, 19.48275862, 19.71929825, 19.92857143, 19.94642857,
           19.96428571, 19.46551724, 19.5       , 19.71929825, 19.73684211,
           20.02777778, 20.        , 19.96428571, 20.        , 19.73684211,
           19.73684211, 19.70175439, 20.01785714, 19.73684211, 19.98214286,
           19.90909091, 19.94642857, 19.73684211, 19.98214286, 19.77192982,
           19.98214286, 19.96428571, 19.44827586, 19.66666667, 19.68421053,
           19.64912281, 19.93877551, 19.94642857, 19.89285714, 19.94642857,
           19.70175439, 19.98214286, 19.91071429, 19.37931034, 19.875     ,
           19.89285714, 19.54054054, 19.75438596, 19.60714286, 19.66666667,
           19.91071429, 19.68421053, 19.92857143, 19.71929825, 19.75438596,
           19.92857143, 19.94545455, 19.68421053, 19.70175439, 19.66666667,
           19.71929825, 19.71929825, 20.01785714, 19.92857143, 19.73684211,
           19.96428571, 19.66      , 19.92857143, 19.43103448, 19.66666667,
           19.94642857, 19.94642857, 19.94642857, 19.96428571, 19.89285714,
           19.71929825, 19.91666667, 19.91071429, 19.92857143, 19.94642857,
           19.64912281, 19.91071429, 19.89285714, 19.89285714, 19.91071429,
           19.71929825, 19.67857143, 19.91071429, 19.92857143, 19.68421053,
           19.68421053, 19.43103448])

  • fp
    (lat, lon, time)
    float32
    0.0 0.0 0.0 ... 2.289e+04 7.629e+03

    array([[[    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
            [    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
            [    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
            ...,
            [    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
            [    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
            [    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ]],
    
           [[    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
            [    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
            [    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
    ...
            [    0.    ,     0.    ,     0.    , ...,  7629.3945,
             22888.184 , 22888.184 ],
            [    0.    ,     0.    ,     0.    , ...,  7629.3945,
              7629.3945, 30517.578 ],
            [    0.    ,     0.    ,     0.    , ...,  7629.3945,
             15258.789 , 22888.184 ]],
    
           [[    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
            [    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
            [    0.    ,     0.    ,     0.    , ...,     0.    ,
                 0.    ,     0.    ],
            ...,
            [    0.    ,     0.    ,     0.    , ...,  7629.3945,
             30517.578 ,     0.    ],
            [    0.    ,     0.    ,     0.    , ...,  7629.3945,
              7629.3945, 22888.184 ],
            [    0.    ,     0.    ,     0.    , ...,     0.    ,
             22888.184 ,  7629.3945]]], dtype=float32)

  • temperature
    (time)
    float32
    15.64 14.35 14.01 ... 13.07 13.31

    units :

    K

    array([15.640076 , 14.350037 , 14.01239  , 14.822693 , 16.491028 ,
           19.916138 , 19.432556 , 17.773743 , 18.512451 , 18.32837  ,
           18.027954 , 16.906433 , 15.878479 , 14.736084 , 14.095093 ,
           13.887695 , 14.175354 , 15.241028 , 17.32959  , 19.320496 ,
           18.878052 , 17.991028 , 17.166199 , 15.750427 , 15.203369 ,
           14.166199 , 13.232178 , 12.398743 , 14.791199 , 16.182556 ,
           17.432556 , 18.6521   , 19.195496 , 19.360779 , 17.878479 ,
           16.144897 , 14.475037 , 13.116211 , 12.712769 , 13.919739 ,
           14.682556 , 16.14447  , 17.896973 , 18.003052 , 18.402344 ,
           18.541138 , 17.15387  , 15.307434 , 14.312317 , 13.475037 ,
           13.147278 , 12.594727 , 13.93811  , 15.075134 , 16.682556 ,
           17.1828   , 15.498413 , 16.453674 , 16.152954 , 15.1796875,
           14.812378 , 15.525574 , 14.333862 , 12.945923 , 13.175354 ,
           15.057556 , 16.243408 , 18.088013 , 16.322449 , 16.44342  ,
           14.60083  , 14.888062 , 14.012695 , 13.27478  , 12.149963 ,
           11.663818 , 12.203857 , 14.699036 , 15.307556 , 15.838013 ,
           17.58551  , 17.463013 , 16.675354 , 14.825012 , 13.850037 ,
           13.621643 , 13.384399 , 13.600037 , 13.816772 , 15.429077 ,
           17.307556 , 19.195496 , 18.610474 , 17.160828 , 16.403687 ,
           15.740662 , 15.300232 , 14.940918 , 15.083862 , 13.078491 ,
    ...
           16.16986  , 15.674622 , 15.980713 , 14.761108 , 16.129272 ,
           15.332764 , 15.860474 , 17.53241  , 18.213013 , 17.81659  ,
           16.703857 , 15.419861 , 15.000427 , 15.050354 , 14.687439 ,
           14.974426 , 14.773743 , 15.291138 , 16.213013 , 17.088013 ,
           18.003052 , 16.97406  , 15.08844  , 15.818298 , 14.241211 ,
           14.269897 , 14.15387  , 13.687439 , 13.356812 , 13.410828 ,
           14.628052 , 15.807556 , 15.901428 , 16.765991 , 16.782837 ,
           14.920532 , 14.2734375, 14.594727 , 15.083984 , 13.777283 ,
           12.393982 , 13.691589 , 17.213013 , 18.122742 , 17.578857 ,
           18.485779 , 18.637451 , 18.802734 , 18.669373 , 17.387695 ,
           16.175354 , 16.041199 , 14.887695 , 15.291138 , 17.253052 ,
           18.445496 , 19.700806 , 20.195496 , 19.894897 , 18.31073  ,
           16.337769 , 15.0824585, 14.312073 , 14.062317 , 15.140381 ,
           15.253052 , 17.100159 , 17.94159  , 15.054749 , 12.9470825,
           13.300354 , 12.753479 , 12.6604   , 12.552673 , 12.800354 ,
           13.007385 , 12.762634 , 13.416321 , 14.695496 , 14.166138 ,
           13.691589 , 12.138611 , 11.544006 , 11.590393 , 11.282776 ,
           11.215393 ,  9.949341 , 10.05658  ,  9.662476 , 11.194763 ,
           12.704041 , 13.350159 , 13.725159 , 14.338013 , 14.075378 ,
           13.074646 , 13.307007 ], dtype=float32)

  • pressure
    (time)
    float32
    994.9 994.3 ... 1.01e+03 1.009e+03

    units :

    hPa

    array([ 994.9192,  994.3181,  994.3651,  994.1363,  993.1855,  992.0133,
            992.0341,  991.4746,  989.7502,  988.8743,  988.4803,  989.3885,
            989.2206,  989.014 ,  989.9035,  990.6757,  991.875 ,  992.3577,
            992.3473,  992.157 ,  992.4304,  992.4638,  992.9948,  993.5303,
            993.6416,  993.4228,  994.0169,  994.2938,  995.1976,  995.1227,
            994.9279,  993.9264,  993.3716,  993.107 ,  993.3469,  993.9938,
            993.6929,  993.4299,  992.8023,  992.5912,  992.5118,  992.4791,
            991.7825,  991.2757,  990.9602,  990.7672,  990.9819,  990.7398,
            990.4157,  990.4048,  990.3318,  990.2883,  990.5333,  990.4781,
            989.9723,  988.5722,  988.2051,  987.6578,  987.0998,  987.4722,
            987.897 ,  987.8259,  987.9325,  987.9838,  989.0248,  989.5853,
            989.918 ,  990.1246,  990.9599,  992.3379,  993.5656,  994.5487,
            995.2004,  996.1456,  996.8638,  998.4511,  999.3455, 1000.374 ,
           1000.97  , 1001.567 , 1001.665 , 1002.171 , 1002.846 , 1003.228 ,
           1003.278 , 1002.925 , 1003.182 , 1003.566 , 1004.327 , 1004.738 ,
           1004.969 , 1004.831 , 1004.648 , 1004.987 , 1005.941 , 1006.241 ,
           1006.819 , 1006.453 , 1006.548 , 1007.446 , 1008.176 , 1008.949 ,
           1008.692 , 1008.611 , 1008.318 , 1008.221 , 1008.988 , 1009.259 ,
           1009.186 , 1008.563 , 1008.86  , 1008.861 , 1009.339 , 1009.201 ,
           1008.925 , 1008.404 , 1007.616 , 1007.736 , 1007.798 , 1007.773 ,
    ...
            997.9325,  997.8781,  997.3772,  997.3143,  997.7148,  997.757 ,
            997.693 ,  997.2031,  997.3762,  997.0297,  996.8853,  996.8566,
            996.3466,  995.451 ,  994.7009,  993.8805,  993.5362,  991.9073,
            990.8405,  990.0166,  988.6914,  988.1125,  988.2752,  988.3303,
            987.9597,  987.1047,  985.6359,  984.8205,  984.6698,  984.4229,
            983.5789,  982.5152,  981.9655,  981.945 ,  982.208 ,  981.8513,
            982.0553,  981.2179,  980.8398,  982.47  ,  982.8874,  984.6092,
            986.1073,  987.1541,  989.037 ,  991.5455,  993.737 ,  995.6596,
            996.8168,  997.8661,  999.0781,  999.5079, 1000.079 , 1000.182 ,
           1000.041 ,  999.3257,  999.0633,  999.1392,  999.002 ,  997.6734,
            996.5509,  995.5034,  994.5828,  992.7145,  992.9497,  991.0938,
            990.456 ,  990.3067,  990.4005,  992.2541,  993.5493,  994.5983,
            995.6412,  995.8322,  996.3949,  996.0748,  996.6905,  996.525 ,
            996.2672,  996.5003,  995.8089,  995.0995,  994.64  ,  994.4482,
            992.8331,  993.1088,  992.7661,  993.1158,  993.8536,  993.8136,
            993.2503,  993.0117,  993.143 ,  993.5527,  994.1947,  994.3773,
            995.4066,  996.7152,  999.4507, 1001.35  , 1003.496 , 1004.382 ,
           1005.46  , 1006.118 , 1006.858 , 1007.723 , 1008.82  , 1009.73  ,
           1009.715 , 1010.123 , 1010.016 , 1009.475 , 1009.695 , 1009.077 ],
          dtype=float32)

  • wind_speed
    (time)
    float32
    8.483 8.844 8.768 ... 3.621 4.743

    units :

    m/s

    array([ 8.483459  ,  8.843811  ,  8.768066  ,  7.3928833 ,  8.15033   ,
            9.448853  ,  9.095276  ,  8.347412  ,  6.9176025 ,  9.319092  ,
            8.4729    ,  8.604614  ,  8.01709   ,  6.5602417 ,  7.1849976 ,
            6.4974365 ,  6.1054077 ,  6.11438   ,  6.8494263 ,  8.208618  ,
           10.104919  ,  7.476074  ,  9.448975  ,  9.109192  ,  8.471924  ,
            8.523376  ,  8.162903  ,  6.4853516 ,  6.623291  ,  7.395691  ,
            6.8895874 ,  6.5928955 ,  8.434448  ,  7.539795  ,  8.404114  ,
            7.553894  ,  8.228821  ,  8.458984  ,  8.828735  ,  7.510498  ,
            5.569275  ,  6.035095  ,  5.5093994 ,  7.680298  ,  5.9473877 ,
            8.0043335 ,  7.420288  ,  7.3426514 ,  8.135376  ,  7.368286  ,
            6.0318604 ,  4.649475  ,  2.1933594 ,  3.0668335 ,  3.5267944 ,
            5.8027954 ,  2.447876  ,  6.1309814 ,  7.4678345 ,  4.0542603 ,
            1.8590088 ,  4.759094  ,  6.744934  ,  5.734375  ,  5.230957  ,
            6.8099365 ,  4.5160522 ,  6.4798584 ,  5.739746  ,  9.626282  ,
            7.1120605 ,  7.6887207 ,  6.8622437 ,  6.33374   ,  6.453491  ,
            7.2006226 ,  7.3204346 ,  8.220825  ,  8.882263  ,  8.351929  ,
            7.678528  ,  7.787964  ,  7.1536255 ,  6.7609863 ,  5.199646  ,
            5.717407  ,  3.7236328 ,  2.3792725 ,  2.8886719 ,  1.0420532 ,
            2.2879028 ,  1.9902344 ,  6.7838745 ,  6.491455  ,  5.1466064 ,
            5.3062744 ,  3.2999878 ,  4.4055176 ,  2.5499878 ,  3.1690674 ,
    ...
            3.1931763 ,  2.4921265 ,  3.140503  ,  5.7384644 ,  4.6742554 ,
            5.30365   ,  6.455139  ,  4.815918  ,  6.2284546 ,  4.5401    ,
            7.494507  ,  8.113525  ,  5.902527  ,  7.2767944 ,  6.9995117 ,
            6.6639404 ,  7.248535  ,  6.059082  ,  6.8433228 ,  6.8933105 ,
            6.8621826 ,  5.408386  ,  7.9352417 ,  6.939514  ,  9.950806  ,
            6.526184  ,  6.4837036 , 11.080078  , 10.697144  ,  9.370056  ,
            8.614258  ,  5.647278  ,  4.967102  ,  1.8279419 ,  1.078125  ,
            6.6557007 ,  6.237671  ,  7.7408447 ,  5.5964966 ,  6.164856  ,
            8.135315  ,  7.5510864 ,  9.898804  , 11.295044  , 10.766663  ,
           10.058472  ,  8.164185  ,  9.777405  , 10.476868  , 10.702759  ,
           10.381165  ,  7.227295  ,  7.5821533 ,  7.7806396 ,  6.902405  ,
            5.586975  ,  7.6796875 ,  8.090393  ,  7.4348145 ,  8.463013  ,
            7.9908447 ,  8.058044  ,  8.801758  ,  9.005676  , 10.714111  ,
            8.63855   , 10.195251  , 10.923157  , 10.746033  ,  9.187561  ,
            9.840332  ,  8.743408  , 10.020325  , 10.070679  , 10.400635  ,
            8.534241  ,  7.5584717 ,  6.960266  ,  4.9536743 ,  5.409485  ,
            8.894165  , 12.123657  , 13.851135  , 12.761414  , 12.367126  ,
           12.106201  ,  9.117859  , 10.151062  , 11.270996  , 10.06311   ,
            9.039673  ,  7.6031494 ,  7.3186646 ,  4.739502  ,  2.4797363 ,
            3.6207275 ,  4.743408  ], dtype=float32)

  • wind_direction
    (time)
    float32
    193.2 203.1 201.4 ... 117.3 172.3

    units :

    degrees

    array([193.18268 , 203.12378 , 201.36737 , 188.14441 , 150.4383  ,
           161.29199 , 178.151   , 170.26758 , 156.27448 , 160.12567 ,
           182.18689 , 208.1311  , 221.38861 , 228.74982 , 236.07098 ,
           231.07391 , 228.8197  , 228.51611 , 201.75763 , 215.32678 ,
           193.0434  , 201.10382 , 213.565   , 203.57928 , 203.1933  ,
           205.05371 , 223.51007 , 222.87158 , 221.71582 , 207.92297 ,
           183.33069 , 186.41193 , 202.4989  , 194.98853 , 206.5415  ,
           204.29071 , 204.40912 , 207.82507 , 199.61963 , 198.71509 ,
           189.52399 , 170.79999 , 170.19861 , 212.8255  , 223.96417 ,
           198.04181 , 191.55951 , 201.1272  , 200.41443 , 197.80573 ,
           215.36682 , 209.43597 , 206.87598 , 161.67639 , 151.7445  ,
           172.4176  , 162.8247  , 137.92212 , 156.67352 , 227.46039 ,
           181.06128 , 217.2052  , 213.17432 , 231.3689  , 233.1756  ,
           258.50403 , 261.45428 , 254.83392 , 277.11243 , 272.68573 ,
           250.37878 , 251.55988 , 244.02051 , 275.79022 , 270.1405  ,
           246.97412 , 247.65472 , 254.1651  , 257.73138 , 271.61603 ,
           252.33118 , 264.44458 , 249.0296  , 247.42072 , 242.90057 ,
           239.75348 , 241.09558 , 255.8222  , 243.87958 , 255.55981 ,
           291.74683 , 296.5011  ,  62.661743,  90.82031 , 118.14862 ,
           126.47223 , 139.97241 , 142.66827 , 155.27728 ,  91.20166 ,
    ...
           188.43292 , 138.8446  , 155.17767 , 145.5904  , 181.29852 ,
           220.10687 , 234.3905  , 219.84271 , 205.55829 , 161.8266  ,
           151.61841 , 179.0567  , 159.90222 , 151.08221 , 171.85638 ,
           185.44019 , 193.91803 , 200.21399 , 216.33917 , 211.42273 ,
           211.7096  , 278.09558 , 223.237   , 248.729   , 249.78918 ,
           298.396   , 294.15552 , 311.8004  , 298.1792  , 297.14362 ,
           311.94177 , 317.95428 , 309.75018 , 335.771   , 112.1322  ,
           107.70032 , 134.91498 , 167.8421  , 210.46777 , 208.21887 ,
           200.0794  , 188.03473 , 177.02173 , 180.90863 , 182.67572 ,
           176.4643  , 174.70349 , 184.09497 , 182.20209 , 201.669   ,
           214.71112 , 222.3028  , 248.11737 , 253.73169 , 262.0052  ,
           234.86438 , 215.69238 , 224.5365  , 202.31073 , 203.03522 ,
           204.69019 , 201.7356  , 197.4909  , 186.76251 , 175.22272 ,
           191.43793 , 205.15527 , 182.22357 , 226.08301 , 207.89948 ,
           207.76031 , 203.25171 , 199.77478 , 196.77551 , 201.18848 ,
           199.48639 , 199.53168 , 220.39581 , 252.49799 , 285.9074  ,
           333.3863  , 343.15143 , 332.27002 , 330.7542  , 329.7298  ,
           334.63342 , 316.78052 , 321.6803  , 322.11603 , 327.64362 ,
           337.5595  , 340.10107 , 337.55798 , 324.05402 ,  42.42151 ,
           117.298706, 172.3429  ], dtype=float32)

  • PBLH
    (time)
    float32
    510.2 341.6 291.5 ... 92.75 92.75

    units :

    m

    array([ 510.18433 ,  341.6247  ,  291.52228 ,  231.89032 ,  510.18433 ,
            686.5085  ,  434.435   ,  475.53497 ,  689.1025  ,  596.3288  ,
            689.1025  ,  357.78168 ,  596.3288  ,  605.9135  ,  476.49268 ,
            510.18433 ,  357.78168 ,  430.66882 ,  519.1636  , 1385.006   ,
           1007.208   , 1264.204   ,  568.7283  ,  543.7831  ,  788.5066  ,
            596.3288  ,  510.18433 ,  420.42578 ,  616.7095  , 1056.7119  ,
            894.5415  , 1081.317   , 1347.792   , 2198.52    ,  510.18433 ,
            788.5066  ,  430.66882 ,  430.66882 ,  482.36072 ,  291.52228 ,
            357.78168 ,  420.7683  , 1007.208   ,  894.5415  , 1039.862   ,
            764.267   , 1057.756   ,  430.66882 ,  545.4724  ,  430.66882 ,
            231.89032 ,  291.52228 ,  312.00012 ,  446.70428 , 1034.931   ,
            510.18433 ,  181.03223 ,  132.50629 ,  894.5415  ,  381.1347  ,
             59.627075,  199.90332 ,  137.29779 ,  178.88507 ,  357.78168 ,
            669.587   ,  785.2084  , 1007.208   ,  132.50629 ,  272.8313  ,
            132.50629 , 1045.431   ,  357.78168 ,  357.78168 ,  291.52228 ,
            329.7068  ,  357.78168 ,  850.8823  ,  510.18433 ,  636.7543  ,
            894.5415  ,  632.9212  ,  438.88373 ,  291.52228 ,  147.36609 ,
            174.09363 ,   92.75372 ,   92.75372 ,  231.89032 ,  614.002   ,
            964.1297  , 1503.621   ,  522.5175  ,  456.86462 ,  357.78168 ,
             88.09845 ,   88.09845 ,   92.75372 ,   59.627075,  103.15802 ,
    ...
             89.33136 ,   59.627075,   59.627075,  167.57593 ,   88.09845 ,
           1169.5449  ,  935.9404  , 1601.185   , 1593.508   , 1734.65    ,
            746.3882  ,  430.66882 ,  618.9514  ,  894.5415  ,  942.2792  ,
           1007.208   ,  689.1025  ,  430.66882 ,  588.1839  ,  768.6522  ,
            731.2208  ,  706.3693  ,  450.0583  ,  240.27051 ,  519.0841  ,
            418.29742 ,  608.3697  ,  943.6366  , 1019.533   , 1009.209   ,
            745.4064  ,  998.9885  ,  658.9596  ,  713.3406  ,  848.5776  ,
            147.36609 ,  894.5415  ,  844.4454  ,  257.87518 ,   92.75372 ,
            178.88507 ,  430.66882 , 1113.996   , 1145.578   ,  774.5368  ,
           1670.046   , 1465.228   , 1126.5071  , 1385.006   ,  689.1025  ,
            942.2792  ,  689.1025  ,  430.66882 ,  497.76642 , 1007.206   ,
            894.5415  , 1524.207   , 1524.208   ,  433.9425  ,  894.5415  ,
            430.66882 ,  450.0583  ,  475.53497 ,  430.66882 ,  231.89032 ,
            510.18433 ,  849.747   , 2138.236   ,  298.36768 ,  720.951   ,
            894.5415  ,  980.5017  ,  894.5415  ,  689.1025  ,  894.5415  ,
            894.5415  ,  596.3288  ,  788.5066  ,  484.70758 ,  596.0555  ,
            979.6506  , 1981.635   , 1926.7639  ,  822.4796  , 1667.28    ,
           1822.521   , 1007.208   ,  872.1522  , 1670.046   , 1097.416   ,
           1285.346   , 1126.5071  , 1755.745   , 1835.145   ,  132.50629 ,
             92.75372 ,   92.75372 ], dtype=float32)

  • release_lon
    (time)
    float32
    1.139 1.139 1.139 ... 1.139 1.139

    units :

    Degrees_east

    array([1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
    ...
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329,
           1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329, 1.1387329],
          dtype=float32)

  • release_lat
    (time)
    float32
    52.52 52.52 52.52 ... 52.52 52.52

    units :

    Degrees_north

    array([52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
    ...
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776,
           52.51776, 52.51776, 52.51776, 52.51776, 52.51776, 52.51776],
          dtype=float32)

  • particle_locations_n
    (height, lon, time)
    float32
    0.0 3.248e-05 0.0 ... 0.0 0.0 0.0

    units :

    long_name :

    Fraction of total particles leaving domain (N side)

    array([[[0.0000000e+00, 3.2484531e-05, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 3.2484531e-05, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 3.2484531e-05, 3.2544136e-05, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            ...,
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             5.1319599e-05, 2.5689602e-05, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 7.7068806e-05, 2.5689602e-05],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             2.5689602e-05, 0.0000000e+00, 0.0000000e+00]],
    
           [[0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [3.2424927e-05, 0.0000000e+00, 3.2544136e-05, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [3.2424927e-05, 0.0000000e+00, 6.5088272e-05, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
    ...
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00]],
    
           [[0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            ...,
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00]]], dtype=float32)

  • particle_locations_e
    (height, lat, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    units :

    long_name :

    Fraction of total particles leaving domain (E side)

    array([[[0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            ...,
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             4.1061640e-04, 3.3390522e-04, 5.3924322e-04],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             3.8498640e-04, 4.3666363e-04, 5.6487322e-04],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             3.5929680e-04, 3.8528442e-04, 3.0815601e-04]],
    
           [[0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
    ...
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00]],
    
           [[0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            ...,
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00]]], dtype=float32)

  • particle_locations_s
    (height, lon, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    units :

    long_name :

    Fraction of total particles leaving domain (S side)

    array([[[0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            ...,
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.]],
    
           [[0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            ...,
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.]],
    
           [[0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            ...,
    ...
            ...,
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.]],
    
           [[0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            ...,
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.]],
    
           [[0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            ...,
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.],
            [0., 0., 0., ..., 0., 0., 0.]]], dtype=float32)

  • particle_locations_w
    (height, lat, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    units :

    long_name :

    Fraction of total particles leaving domain (W side)

    array([[[0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            ...,
            [3.2424927e-05, 9.7453594e-05, 6.5088272e-05, ...,
             0.0000000e+00, 2.5689602e-05, 0.0000000e+00],
            [3.2424927e-05, 3.2484531e-05, 0.0000000e+00, ...,
             2.5689602e-05, 0.0000000e+00, 0.0000000e+00],
            [6.4909458e-05, 0.0000000e+00, 3.2544136e-05, ...,
             2.5689602e-05, 0.0000000e+00, 0.0000000e+00]],
    
           [[0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
    ...
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00]],
    
           [[0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            ...,
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00],
            [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,
             0.0000000e+00, 0.0000000e+00, 0.0000000e+00]]], dtype=float32)

  • mf_mod
    (time)
    float32
    37.35 35.34 27.8 ... 14.56 24.52

    resample_to :

    coarsest

    array([37.34853  , 35.33773  , 27.799513 , 13.802546 ,  4.2526965,
            2.971444 ,  3.8846684,  3.845897 ,  3.5672452,  8.569703 ,
           18.37524  , 25.673927 , 23.749813 , 29.048111 , 36.23223  ,
           29.951513 , 30.297514 , 21.976889 , 12.020665 ,  9.728594 ,
            9.67127  , 13.30284  , 19.33358  , 25.79167  , 27.717373 ,
           25.99897  , 27.062805 , 21.176718 , 17.307856 , 12.722499 ,
            6.865408 ,  8.830698 , 11.165742 , 12.265608 , 19.983404 ,
           26.795385 , 29.328102 , 31.21886  , 27.34428  , 14.585773 ,
           12.061527 , 12.898601 , 12.370915 , 11.179406 , 11.780992 ,
           14.2984915, 20.55703  , 26.487282 , 31.10362  , 34.27147  ,
           37.182205 , 35.444717 , 29.228277 , 25.327671 , 21.343357 ,
           15.200205 , 12.881175 , 11.643733 , 19.447935 , 33.918385 ,
           43.415062 , 49.443954 , 56.77151  , 54.422287 , 53.415028 ,
           46.063488 , 38.12833  , 35.646805 , 33.95901  , 36.16345  ,
           40.796173 , 45.62204  , 47.462433 , 43.67146  , 39.41717  ,
           22.607574 , 15.732253 , 19.869953 , 19.210604 , 14.738877 ,
           14.989182 , 17.362026 , 19.427746 , 23.6137   , 26.174051 ,
           29.806831 , 36.14797  , 44.322674 , 37.430424 , 32.525494 ,
           26.527906 , 25.975939 , 28.052998 , 27.961708 , 22.275225 ,
           20.369507 , 24.274574 , 31.64384  , 36.470432 , 36.308113 ,
    ...
           21.407494 , 19.289438 , 20.201641 , 28.795406 , 45.780743 ,
           44.194874 , 21.895784 , 12.528505 , 15.709129 , 12.5210495,
            7.648608 ,  8.7779045, 10.301787 , 11.912912 , 16.128124 ,
           27.956572 , 29.975779 , 27.17092  , 22.508295 , 19.518105 ,
           19.265263 , 29.20205  , 37.95838  , 42.959415 , 63.64801  ,
           61.85968  , 19.619587 , 16.562286 , 13.062411 , 13.529734 ,
           10.382277 ,  6.265456 ,  5.5641775,  7.375521 ,  9.886796 ,
            8.831683 , 14.695633 , 37.345837 , 41.560703 , 47.358315 ,
           31.161825 , 10.963774 ,  5.781119 ,  5.608469 ,  7.3720646,
            8.916411 , 11.848292 , 14.561811 , 15.291674 , 21.16     ,
           24.163527 , 26.593231 , 16.355646 , 13.170094 , 10.508017 ,
            9.208046 , 10.592859 , 14.017923 , 16.56197  , 20.370354 ,
           25.675346 , 30.14644  , 31.669052 , 22.811531 ,  9.588773 ,
            7.924422 ,  8.085369 ,  9.32877  , 10.092244 , 11.178455 ,
           12.982593 , 18.364685 , 24.049093 , 27.205849 , 27.529284 ,
           28.224747 , 23.952017 , 25.349268 , 30.865599 , 35.693157 ,
           18.691717 ,  3.0909324,  3.4081225,  3.6038013,  4.254746 ,
            5.1437197,  5.6772933,  6.2337785,  3.1979098,  3.488845 ,
            4.784305 ,  5.8790836,  5.207522 ,  7.2529335, 12.187598 ,
           14.557251 , 24.51839  ], dtype=float32)

• Attributes: (40)

  author :

  OpenGHG Cloud

  processed :

  2022-03-22 10:42:55.132208+00:00

  data_type :

  footprints

  site :

  tac

  height :

  100m

  domain :

  europe

  model :

  NAME

  start_date :

  2012-08-01 00:00:00+00:00

  end_date :

  2012-08-31 22:00:00+00:00

  max_longitude :

  39.38

  min_longitude :

  -97.9

  max_latitude :

  79.057

  min_latitude :

  10.729

  time_resolution :

  standard_time_resolution

  heights :

  [ 500. 1500. 2500. 3500. 4500. 5500. 6500. 7500. 8500. 9500. 10500. 11500. 12500. 13500. 14500. 15500. 16500. 17500. 18500. 19500.]

  variables :

  ['fp', 'temperature', 'pressure', 'wind_speed', 'wind_direction', 'PBLH', 'release_lon', 'release_lat', 'particle_locations_n', 'particle_locations_e', 'particle_locations_s', 'particle_locations_w']

  data_owner :

  Simon O'Doherty

  data_owner_email :

  s.odoherty@bristol.ac.uk

  inlet_height_magl :

  100m

  comment :

  Cavity ring-down measurements. Output from GCWerks

  long_name :

  tacolneston

  instrument :

  picarro

  sampling_period :

  60

  inlet :

  100m

  port :

  9

  type :

  air

  network :

  decc

  species :

  ch4

  conditions_of_use :

  Ensure that you contact the data owner at the outset of your project.

  source :

  In situ measurements of air

  Conventions :

  CF-1.8

  file_created :

  2022-03-22 10:42:54.957441+00:00

  processed_by :

  OpenGHG_Cloud

  sampling_period_unit :

  s

  station_longitude :

  1.13872

  station_latitude :

  52.51775

  station_long_name :

  Tacolneston Tower, UK

  station_height_masl :

  50.0

  scale :

  WMO-X2004A

  resample_to :

  coarsest

When the same calculation is being performed for multiple methods, the last calculation is cached to allow the outputs to be produced more efficiently. This can be disabled for large datasets by using cache=False.

For a ModelScenario object, different analyses can be performed on this linked data. For example if a daily average for the modelled observations was required, we could calculate this setting our resample_to input to "1D" (matching available pandas time aliases):

modelled_observations_12H = scenario.calc_modelled_obs(resample_to="1D")
modelled_observations_12H.plot()

Caching calculated data

[<matplotlib.lines.Line2D at 0x7f0990631910>]

To allow comparisons with multiple fluxes inputs, more than one flux source can be linked to your ModelScenario. This can be either be done upon creation or can be added using the add_flux() method. When calculating modelled observations, these flux sources will be aligned in time and stacked to create a total output:

scenario.add_flux(species=species, domain=domain, source=source_energyprod)

Adding flux to model scenario

scenario.plot_comparison()

Caching calculated data

Output for individual sources can also be created by specifying the sources as an input:

# Included recalculate option to ensure this is updated from cached data.
modelled_obs_energyprod = scenario.calc_modelled_obs(sources="energyprod", recalculate=True)
modelled_obs_energyprod.plot()

Caching calculated data

[<matplotlib.lines.Line2D at 0x7f0990634fd0>]

Plotting functions to be added for 2D / 3D data

---

previous

Workflow 2: handling overlapping observation data

next

Workflow 4: Working with carbon dioxide